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ABSTRACT 


Pease  (1987)  related  the  effects  of  atmospheric  forcing,  mainly  temperature  and  wind 
fields,  to  the  size  of  coastal  polynas.  Using  Pease's  formulation  and  Kilhvorth's  (1977) 
plume  model  as  applied  by  Melling  and  Lewis  (19S2),  the  effects  of  atmosoheric  forcing 
on  brine  injection  into  the  Arctic  pycnoeline  are  investigated.  This  paper  will  discuss  the 
likelihood  of  coastal  polynas  as  a  source  for -denser  abyssal  waters. 

A  standard  case  was  developed  for  the  model  with  initial  conditions  taken  from 
Melling  and  Lewis  (19S2)  and  Pease  (1987)  for  comparison  with  individual  sensitivity 
experiments.  Ten  environmental  parameters  were  individually  examined  for  their  influ¬ 
ence  on  the  plume  depth  after  90  days.  The  standard  case  resulted  in  a  90-day  plume 
depth  of  436  meters.  A  submarine  canyon  case  was  simulated,  resulting  in  plume  pene¬ 
tration  to  over  1300  meters  in  90  days.  Further  experiments  used  actual  T-S  soundings 
from  Aagaard  et  al.  (19S1)  and  Ostlund  et  al.  (1987).  Finally,  a  20  kilometer  wide  plume 
is  shown  to  penetrate  to  almost  600  meters  in  90  days. 
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I.  INTRODUCTION 


The  Arctic  Ocean,  a  small  and  remotely  located  part  of  the  world's  ocean,  has  de¬ 
manded  the  attention  of  those  who  wish  to  understand  the  Earth's  climate.  The  renewed 
emphasis  in  recent  years  given  to  understanding  the  Earth's  climate  has  underscored  the 
role  of  the  Arctic  Ocean  in  the  scheme  of  the  global  energy  balance.  The  Arctic  Ocean 
is  also  a  key  player  in  the  world  oceans'  thermohaline  circulation.  Because  of  the  im¬ 
portance  of  the  Arctic  Ocean  in  the  global  picture,  understanding  the  nuances  of  the 
Arctic  Ocean's  physical  oceanography  and  its  sensitivities  to  environmental  change  is 
essential  before  climatologists,  oceanographers,  and  meteorologists  can  develop  a  more 
realistic  and  longer  range  global  climatological  concept.  Two  puzzling  aspects  of  the 
Arctic  Ocean  are  its  isothermal  halocline  and  the  more  saline  deep  waters.  How  is  the 
halocline  maintained  in  a  nearly  horizontally  independent  state?  What  is  the  mechanism 
behind  the  formation  of  the  Arctic  Ocean's  deep  waters?  These  topics  will  be  discussed 
in  this  thesis. 

Throughout  the  Arctic  Ocean,  the  characteristics  of  the  upper  layers  are  remarkably 
independent  of  the  horizontal  location  of  the  sample.  The  water  column  consists  of  a 
mixed  layer,  halocline,  thermocline,  and  deeper  waters.  The  mixed  layer  is  a  relatively 
fresh  water  layer  approximately  30  meters  deep.  River  run-off  and  the  summer  cycle  of 
melting  ice  are  considered  the  main  contributors  to  the  fresh  water  layer.  The  halocline 
is  the  layer  of  water  beneath  the  mixed  layer.  In  this  layer,  salinities  rise  from  31  to  as 
high  as  34.6  ppt  at  the  200  meter  depth  level  ( Kill  worth  and  Smith,  1984).  The  halocline 
is.  in  essence,  an  isothermal  layer  with  temperatures  just  above  the  freezing  point.  The 
thermocline  is  a  layer  below  the  halocline.  in  which  temperatures  rise  to  approximately 
0°  C  at  250  meters  (Melling  and  Lewis.  1982)  and  continues  to  rise  to  0.5s  C  at  400  meters 
due  to  the  layer  of  warm  Atlantic  water.  Below  this  depth,  temperatures  decrease  with 
increasing  depth.  The  deep  waters  of  the  Arctic  Ocean  are  among  the  most  dense  waters 
found  in  the  world's  oceans  and  have  been  suggested  as  being  an  additional  source  for 
ventilating  the  world's  oceans'  deep  thermohaline  circulation  (Aagaard,1981).  Although 
the  physical  oceanography  of  the  Arctic  Ocean  is  well  described  in  literature,  the  exact 
mechanism(s)  for  the  maintenance  of  the  halocline  and  deep  water  formation  are  not 
fully  understood. 
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The  halocline  is  generally  isothermal  except  in  the  Canadian  Basin,  where  warm 
Pacific  waters  entering  via  the  Bering  Strait  causes  a  warm  layer  at  about  the  75-100 
meter  depths,  and  presents  a  paradox  when  examining  its  origin.  Sandwiched  between 
a  cold,  relatively  fresh  layer  and  a  warm,  saline  layer,  an  obvious  solution  to  the  main¬ 
tenance  of  the  halocline  might  be  simple  vertical  mixing  between  the  two  layers.  How¬ 
ever,  vertical  mixing  of  the  cold,  less  saline  surface  water  with  the  warmer,  more  saline 
Atlantic  water  is  not  considered  to  be  a  source  since  heat  is  not  mixed  along  with  salt. 
Aagaard  et  al.  (1981)  showed  a  discontinuity  existed  in  the  T-S  correlation  slope  at 
about  150  meters.  Either  salinization  of  surface  waters  by  freezing  (i.e.  brine  rejection) 
or  cooling  of  Atlantic  water  could  be  two  possible  mechanisms  to  produce  halocline 
water.  The  authors  (Aagaard.  et  al.,  1 9 S 1 )  demonstrate  mixing  between  the  two  layers 
cannot  be  a  source  of  halocline  water  by  examining  the  required  salinity  of  the  surface 
water.  Using  a  T-S  diagram,  and  assuming  linear  mixing  at  a  ratio  of  3:2  (surface  water 
to  Atlantic  water),  the  surface  water  would  need  to  have  a  salinity  of  34.7  ppt.  Salinity 
this  high  has  not  been  observed  in  the  Arctic  surface  layer.  Therefore,  the  maintenance 
of  the  Arctic  halocline  is  generally  thought  to  be  attributable  to  lateral  advection  of  cold, 
saline  water  from  the  continental  shelves.  The  shelf  waters  are  made  more  saline  due  to 
brine  rejection  associated  with  ice  growth.  Aagaard  (1981)  proposed  that  continental 
shelf  areas  where  an  annual  growth  of  ice  less  than  2  meters  is  necessary  to  modify 
summertime  salinity  profiles  to  wintertime  profiles  may  be  the  source  areas  of  the 
halocline  waters.  Upwelling  of  Atlantic  water  onto  the  shelves  and  subsequent  cooling 
has  also  been  suggested  as  a  possible  source.  Upwelling  along  submarine  canyons  has 
been  reported  as  early  as  1929  (Sverdrup,  1929).  As  further  evidence  eluding  to  the 
halocline's  possible  shelfwater  origin,  water  on  the  continental  shelves  has  been  shown 
to  have  temperature  and  salinity  characteristics  similar  to  the  waters  found  offshore  at 
a  depth  of  150  meters  (Melling  and  Lewis,  1982). 

In  their  paper,  Melling  and  Lewis  (19S2)  discuss  the  lateral  advection  of  low  tem¬ 
perature  high  salinity  shelfwaters  into  the  halocline,  with  particular  attention  paid  to  the 
effects  of  brine  rejection  due  to  the  growth  of  ice.  Aagaard  (1981)  used  two  methods 
of  computing  the  necessary  production  rate  of  high  salinity  shelfwater  to  maintain  the 
halocline.  By  considering  the  production  rate  required  to  renew  the  halocline  in  10  years 
and  necessary  mixing  ratios  to  achieve  desired  T-S  characteristics,  Aagaard  suggests  a 
2.5.y106  to  4.0.vl0-m3  —  sec-1  for  a  production  rate  of  cold  saline  shelfwater,  assuming  a 
halocline  renewal  period  of  10  years.  In  order  to  increase  the  ambient  shelfwater  salinity 


to  values  which  would  form  sufficiently  dense  water  to  descend  the  slope  into  the 
halocline,  Melling  and  Lewis  (1982)  emphasize  the  important  role  of  ice  dynamics,  spe¬ 
cifically  the  role  of  divergence.  Areas  of  large  wintertime  ice  cover  divergence  allow 
rapid  ice  growth  hence  increased  salinization  of  the  underlying  water  column.  Leads  and 
polynas  are  locations  of  rapid  ice  growth  due  to  the  divergence  of  the  insulating  ice 
cover. 

Schumacher  et  al.  (  1983)  examined  the  degree  of  salinization  of  the  underlying  water 
column  beneath  a  polyna.  In  their  study,  measurements  of  salinity  and  water  current 
velocity  beneath  a  polyna  located  at  St.  Lawrence  Island  indicated  polynas  to  be  an  area 
of  high  ice  growth  rates  resulting  in  density  driven  currents  due  to  increased  salinization 
(Schumacher,  et  al..  1983).  Coupling  the  salt  source,  i.c.  polyna,  with  the  interleaving 
plumes  of  dense  shelfwater  found  by  Melling  and  Lewis  (1982)  may  provide  a  possible 
mechanism  for  the  maintenance  of  the  Arctic  halocline.  The  polyna  provides  an  atmo¬ 
spheric  dependent  source  of  salt  into  the  water  column  possibly  forming  plumes  which 
descend  down  into  the  halocline,  as  described  by  Melling  and  Lewis  (19S2).  Along  with 
the  mechanism  associated  with  the  halocline's  maintenance,  the  formation  of  the 
Canadian  Basin's  deep  water  has  also  been  attributed  to  having  a  source  of  low  tem¬ 
perature,  high  salinity  shelfwater. 

The  dense  waters  which  ventilate  the  world’s  oceans  flow  out  of  the  Arctic  Ocean 
over  the  Greenland-Scotland  ridge  system  at  a  depth  of  about  2600  meters.  A  single 
source  was  advocated  by  Helland-Hansen  and  Nansen  (1909)  in  which  deep  water  is 
formed  in  the  Greenland  Sea.  Nansen  ( 1906)  had  written  earlier  of  the  possibility  of  the 
Greenland  Sea  Deep  Water  being  modified  by  cooling  in  the  Barents  Sea.  This  theory 
omitted  the  possibility  of  deep  water  formation  within  the  two  basins  of  the  Arctic 
Ocean,  namely  the  Eurasian  and  Canadian  Basins.  In  1981,  Aagaard  postulated  that 
deep  water  is  also  formed  within  the  Arctic  Ocean  due  to  the  salinity  distribution. 
Greenland  Sea  Deep  Water  has  a  salinity  of  approximately  34.90  ppt.  which  is  the  least 
saline  of  Arctic  deep  waters.  Therefore,  at  least  two  sources  of  deep  water  formation 
are  thought  to  exist.  One  source  is  the  traditional  Greenland  Sea  formation  area  and  the 
new  source  proposed  by  Aagaard,  is  within  the  Arctic  Ocean's  two  basins. 

The  characteristics  of  water  mass  properties  found  in  the  Canadian  and  Eurasian 
Basins  are  similar  until  about  the  depth  of  the  Lomonosov  Ridge.  At  about  1500  meters 
depth,  the  Canadian  Basin  becomes  more  saline  and  is  slightly  warmer  than  the  Eurasian 
Basin  Deep  Water  (Aagaard.  1981).  Under  the  traditional  viewpoint  of  Greenland  Sea 
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Deep  Water  Tilling  the  two  basins  through  lateral  advection,  salinities  should  not  exceed 
34.9  ppt.  However,  during  the  LOREX  (Lomonosov  Ridge  [Experiment;  Weber,  1979) 
salinities  in  the  b  .:,ns  vvere  determined  to  be  greater  than  34.9  ppt.  This  suggests  an 
additional  source  of  salt  if  not  a  completely  new  source  of  deep  water.  Current  meters, 
placed  o' \'r  the  Lomonosov  Ridge,  gave  evidence  to  the  existence  of  some  tidal  and 
other  episodic  surges  carrying  salini/cd  Greenland  Sea  Deep  Water  over  the  ridge  into 
the  Canadian  Basin.  Even  with  the  surges  over  the  Lomonosov  Ridge  into  the  Canadian 
Basin,  water  from  the  Eurasian  Basin  must  have  an  additional  salt  source  to  obtain  the 
characteristics  of  water  found  in  the  Canadian  Basin  (Aagaard,  1981).  Residence  times 
can  be  calculated  by  comparing  the  mixing  ratios  of  shelf  water  to  Atlantic  Intermediate 
water  necessary  to  achieve  the  T-S  characteristics  found  in  the  Canadian  Basin  Deep 
Water.  Considering  a  small  flux  of  shelfwater  (.003  to  .06  Sv)  necessary  to  ventilate  the 
Canadian  Basin  over  several  hundred  years  may  also  explain  the  warmer  temperatures 
found  in  the  Canadian  Basin  (Aagaard  ct  ah,  1985).  By  using  tracers  such  as  (5lf,0  and 
l4C,  Ostlund  et  al.  (1987)  agree  with  the  magnitudes  of  the  renewel  time  scales  for  both 
basins  presented  by  Aagaard  et  al.(  1 985).  Ostlund  et  al.  (1987)  and  Aagaard  et  al.  (1985) 
place  the  residence  time  of  the  Eurasian  Basin  in  the  order  of  decades  vice  several  hun¬ 
dred  years  (700  years,  Ostlund,  ct  ah,  1987)  for  the  Canadian  Basin.  The  requirement 
remains  of  an  additional  salt  source  at  small  fluxes  to  ventilate  the  Canadian  Basin  over 
several  hundred  years  point  to  an  episodic  mechanism.  Killworth  and  Smith  (1984)  write 
that  shelfwater  sinking  at  the  basin  edges  may  occur  but  a  source  of  waters  sufficiently 
dense  to  penetrate  the  Arctic  Ocean's  stable  mid-layer  remains  a  problem.  Again,  using 
a  coastal  polyna  as  a  salt  source,  could  shelfwater  be  modified  under  extreme  conditions 
to  produce  such  dense  water  which  would  then  descend  to  the  deeper  depths  via  gravity 
flow  plumes? 

This  thesis  will  combine  Pease's  numerical  formulations  relating  atmospheric  forc¬ 
ing,  mainly  wind  and  air  temperature,  to  the  opening  and  closing  of  coastal  polynas,  with 
a  streamtube  gravity  flow  model  developed  by  Smith  (1975).  With  simplified 
hydrography  and  suitable  environmental  conditions,  this  model  can  produce  events  of 
deep  plume  penetration  along  with  a  more  consistent  source  of  high  salinity,  low 
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temperature  water  for  the  maintenance  of  the  Arctic  Ocean's  halocline.  Chapter  2  will 
describe  the  numerical  model  used  in  this  study  along  with  an  examination  into 
drawbacks  associated  with  Smith's  streamtube  model  as  described  by  other  studies,  i.e. 
Smith  (1975).  Killworth  (1977).  and  Melling  and  Lewis  (19S2).  Chapter  3  will  describe 
results  obtained  by  altering  environmental  inputs  into  the  model  and  Chapter  4  will 
discuss  the  results  and  implications. 


II.  NUMERICAL  METHODS  AND  PARAMETERS 


A.  THE  NUMERICAL  MODEL 

With  the  increased  desire  to  understand  the  nature  of  the  Arctic  halocline  and  its 
origin,  there  existed  sufficient  pieces  of  the  puzzle  to  develop  a  coupled  model.  The  3-D 
turbulent  entraining  plume  model  developed  by  Smith  (1975)  has  been  used  by  several 
studies  (Killworth  ( 1977).  Kilhvorth  and  Carmack  ( 197S).  and  Melling  and  Lewis  ( 1982)) 
to  gain  insight  into  gravity  flow  currents  arising  from  brine  rejection  into  rhe  water  col¬ 
umn  due  to  ice  growth.  Melling  and  Lewis  (19S2),  compute  the  initial  density  of  the 
plume  by  estimating  the  average  ice  grown  annually  in  a  shallow  shelf  and  and  the  re¬ 
sultant  salinity  contribution.  In  this  study.  Pease's  (19S7)  formulation  relating  a 
polyna's  size  to  offshore  wind  speed  and  air  temperature,  provides  a  specific  source  of 
brine  which  is  then  coupled  with  Smith's  plume  model,  as  used  by  Melling  and  Lewis 
( 19S2),  in  order  to  study  the  possible  contribution  of  the  polyna  produced  denser  waters 
to  the  Arctic  halocline  and  possible  formation  of  Canadian  Basin  Deep  Water. 

1.  The  coastal  polyna 

Pease  (1987)  describes  the  polyna's  maximum  width  as  an  equilibrium  between 
frazil  ice  production  and  its  advection  rate.  The  frazil  ice  production  rate  is  a  function 
of  the  heat  exchange  at  the  air-sea  interface.  Therefore,  with  wind  speed  and  air  tem¬ 
perature  as  independent  variables,  one  can  compute  the  size  of  the  polyna  and  its  frazil 
ice  production  rate. 

The  polyna  size  can  be  found  by  solving: 


(IX 

dt 


=  Vi~ 


XpF, 

II; 


(1) 


where  is  the  change  in  polyna  width  with  time.  V,  is  the  advection  rate  of  frazil  ice, 
X.  is  the  polyna  width,  F,  is  the  frazil  ice  production  rate,  and  H,  is  the  collection  depth 
of  grease  ice  or  nilas.  The  solution  of  equation  (1)  assumes  the  meteorological  condi¬ 
tions  are  steady  throughout  the  polyna's  opening  time,  thus  allowing  treatment  as  a 
linear  differential  equation.  Furthermore,  if  the  meteorological  event  is  assumed  to  be 
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of  sufficient  duration  so  that  the  polyna  reaches  maximum  width,  equation  (1)  can  be 
written  as: 


,  I  11; 

X .{maximum )  =  — j—  .  (2) 

*  i 

This  allows  solving  equation  (2)  fox  the  maximum  width  after  computing  the 
frazil  ice  production  rate.  Fr  In  this  study,  meteorological  events  of  interest  are  4-10 
days  duration,  therefore  this  assumption  allowing  equation  (2)  is  valid.  Ou  (19SS)  de¬ 
scribes  a  more  detailed  coastal  polyna  model  and  concludes  the  ice  edge  is  less  affected 
by  higher  frequency  atmospheric  disturbances  than  longer  period,  synoptic  type  vari¬ 
ations.  Although  Ou  (1988)  includes  additional  physics  to  model  the  temporal  changes 
of  the  polyna's  ice  edge,  the  steady  state  solution  is  still  that  as  given  by  Pease  (1987) 
and  shown  above  in  equation  (2). 

The  frazil  ice  production  rate  is  given  by  Pease  (1987)  as: 


lr7egTg  ~  Qlu  +  Pa£h£pK^  Tg  ~ 

PiL 


(3) 


where  the  evaporative  heat  flux  has  been  neglected  due  to  its  small  contribution  relative 
to  the  uncertainty  of  the  sensible  heat  flux.  Q,u  is  the  upward  longwave  radiation,  aej^ 
is  the  downward  longwave  radiation,  L  is  the  latent  heat  of  freezing  for  salt  water,  p,  is 
the  density  of  young  sea  ice.  and  ppChCpV£Ta  —  T„)  is  the  sensible  heat  flux. 

2.  Salinity  distribution  within  the  water  column 

The  salinity  contribution  to  the  water  column  from  frazil  ice  growth  is  assumed 
to  be  uniform  with  depth  and  is  computed  using: 


—  (f>  —  St)F,A. 

S  =  S'  + - ,  —  1  ' 

^new  °  < 


{  qdz) 


(4) 
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from  Killworth  and  Smith  (19S4).  Srew  is  the  new  salinity  value  of  the  water  column  after 
brine  rejection.  5  is  the  average  salinity  of  the  shelfwater,  5.  is  the  salinity  of  the  ice,  A. 
is  the  effective  area  of  salinization,  and  f  qdz  is  the  areal  flux  of  the  underlying  water. 
Using  F,  from  the  polyna  model,  one  can  compute  the  salinity  of  shelfwater  exiting  a 
polyna. 

3.  Equation  of  state 

The  density  of  the  current  is  computed  using  the  International  one  atmosphere 
equation  of  state  developed  by  Millero  and  Poisson  (1981).  The  equation  is  given  below. 

p  =  p0  4-  AS  4-  BS  4-  CS  (5) 


where, 

A  =  8.24493  x  10"’  -  4.0899  x  10_3r  4-  7.6438  x  10_5r2  -  8.2467  x  10"  V  +  5.3875  x  10'V 
B  =  -5.72466  x  10“3  4-  1.0227  x  10~4r  -  1.6546  x  10  -  6r2 
C  =  4.83 14  x  10'4 

and  p0  is  the  standard  density  of  seawater.  This  allows  the  calculation  of  the  current's 
relative  buoyancy  using: 


A. 


A  Po  ’ 


(9) 


In  this  fashion,  the  results  from  the  polyna  model  can  be  used  as  inputs  to  the 
streamtube  model. 

4.  The  streamtube  model 

Smith's  (1975)  model  of  boundary  currents  includes  the  effects  of  friction, 
entrainment,  and  Coriolis  forces.  The  3-D  plume  was  modified  by  Killworth  (1977)  by 
altering  the  entrainment  rate.  Killworth  (1977)  related  the  entrainment  rate  to  the 
product  of  the  plume's  area  and  velocity  vice  only  the  velocity  (Smith,  1975).  Melling 
and  Lewis  ( 19S2)  used  Killworth's  version  of  the  plume  model  with  an  entrainment  co¬ 
efficient  dependent  on  the  Richardson  number.  The  model  used  in  this  paper  is  the 
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Melling  and  Lewis  (1982)  adaptation  with  the  inclusion  of  several  stratified  layers  to  al¬ 
low  a  single  model  initialization. 

Smith's  plume  model  uses  two  coordinate  systems.  The  Cartesian  system  is 
aligned  with  the  bottom  so  the  X-axis  lies  along  the  shore  and  the  Y-axis  lies  along  the 
slope.  The  second  coordinate  system  is  the  curvilinear  system  with  2  and  ;/  axes.  The 
position  of  the  plume  axis  can  be  defined  by  the  value  of  c  .  The  angle  which  the  plume 
crosses  the  isobaths  is  given  by  the  angle  /?  .  The  coordinate  system  is  illustrated  in 
Figure  1  on  page  10  and  includes  gravitational,  g,  and  rotational,  Q,  vectors. 

The  equations  for  the  path  flow  aref 


(10) 


and 


(11) 


The  plume  is  considered  to  be  a  well-mixed,  broad,  thin  layer  of  high  density  fluid  adja¬ 


cent  to  the  bottom.  Pressure  gradient  forces  pull  the  plume  down  the  slope  while 
Coriolis  forces  bend  the  flow  to  the  right.  If  frictional  forces  were  not  considered,  the 
flow  would  move  horizontally  along  the  isobaths  (Killworth.  1977)  in  geostrophic  bal¬ 
ance.  Friction  is  included,  using  quadratic  drag  laws  to  relate  frictional  resistance  to  the 
square  of  the  mean  velocity  (Smith,  1975).  The  plume  is  thus  allowed  to  flow  down  the 
slope,  constantly  entraining  ambient  fluid  and  losing  its  negative  buoyancy.  Once  the 
plume  reaches  neutral  buoyancy,  the  plume  interleaves  with  the  ambient  waters.  The 
equations  of  the  plume's  flow  from  Killworth  (1977),  are: 


^r{AV)  =  E,A^V, 

“S 


(13) 


-jT-  (A  VA)  =  -  A  VN1  sin  Q  sin  /?, 


(14) 
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A  sin  9  cos  P  —  JV. 


(15) 


These  equations  are,  the  along-stream  derivatives  of  the  mass  flux  (12),  the  momentum 
flux  (13),  the  buoyancy  flux  (14)  and  the  cross  stream  momentum  balance  (15).  In  the 
above  equations,  A  is  the  cross-sectional  area  of  the  plume,  V  is  the  plume's  velocity,  A 
is  the  plume  s  buoyancy,  9  is  the  slope  angle,  p  is  the  angle  at  which  the  plume  crosses 
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the  isobath,  N  is  the  Brunt- Viiisala  frequency,  Kx  is  the  drag  coefficient, /is  the  Coriolis 
parameter,  and  £,  is  the  entrainment  coefficient. 

The  equations  were  integrated  using  a  modified  Adams  method. 

Predictor: 


Corrector: 


37+ 1  =  >7  +  "Tf  (^/  -  59/_,  +  37 'fi_2  -  9/_3) 


(16) 


37+i  =>V  +  “ 24  (9/+1  +  19/-  o/_,  +/_2), 


(16a) 


with  the  initial  steps  provided  with  the  Runge-Kutta  method, 

>7+i  =>'i  +  y  (*t  +4/S:2  +^3) 

*1  =A-X ,■•})) 

^2  =Axi  +  y  •  >7  +  y  ^1) 

^3  =/(*,- + ^->v ~ m») 


(IS) 


This  model  uses  Killworth's  (1977)  entrainment  and  frictional 

parameterizations,  which  has  the  entrainment  rate  and  frictional  force  porportional  to 
the  product  of  the  area  and  velocity  (velocity  squared  for  friction).  Thus, 


E 

=  ( —  )t ,  where  w  is  the  width  of  the  plume  and  b  is  the  depth  of  the  plume. 
E\ 


K 


w  ,  1 


Also  used  in  this  model  is  Bo  Pedersen's  (19S0)  entrainment  parameterization  for  small 
slopes.  Melling  and  Lewis  (19S2)  included  this  in  their  version  of  Smith's  model.  The 

E  =  0.072  sin  6  sin  p  =  0.072)  -§r )  where 

Ki 


entrainment  coefficient  is  given  as 
Ri  =  bA  cos  d/V2  is  the  bulk  Richardson  number  of  the  How. 

5.  The  coupled  model 

The  model  takes  the  form  illustrated  in  Figure  2  on  page  13.  Underneath  the 
polyna,  a  plume  with  a  depth,  b.  and  a  width,  w,  is  salinized  due  to  brine  rejection.  The 
plume  has  an  initial  velocity,  V„  .  The  effects  of  varying  the  initial  velocity  will  be  ex¬ 
amined  in  Section  3,  however,  the  standard  case  value  of  V„  =  .OAm/s  is  used.  This  value 
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is  dose  to  the  values  used  by  Melling  and  Lewis  (19S2)  and  observed  values  of  currents 
adjacent  to  polynas  reported  by  Schumacher  et  al.  (19S3).  The  value  of  V0  determines 
the  areal  flux  used  in  calculating  the  salinity  increase  in  the  water  column.  The  initial 
velocity  is  increased  after  the  non-entraining,  non-rotating  rudiment  plume  travels  10 
km  straight  downslope.  Effects  of  rotation  are  neglected  at  this  early  stage  of  the  plume 
to  minimize  the  exposure  to  brine  rejection  under  the  polyna.  Increasing  the  velocity 
after  10  km  ensures  more  realistic  salinity  values  of  the  plume,  and  would,  in  any  case, 
occur  in  the  first  steps  of  the  streamtube  model.  The  non-entraining,  non-rotating 
rudiment  plume  continues  straight  downslope  until  a  depth  of  45  meters  is  reached.  This 
depth  was  chosen  due  to  its  use  by  Melling  and  Lewis  (19S2)  as  the  initial  depth  of  the 
streamtube  model. 

At  45  meters,  the  plume  descends  the  slope  according  to  the  dynamics  set  forth 
and  established  by  Smith  (1975),  Killworth  (1977).  and  Melling  and  Lewis  (19S2).  Three 
stratified  layers  are  defined  in  this  coupled  model.  Layer  1  lies  between  45  and  55  meters 
depth.  Layer  2  lies  between  55  and  300  meters,  and  Layer  3  lies  below  300  meters.  The 
plume  is  assumed  to  interleave  with  ambient  waters  when  the  velocity  is  less  than  .01 
m's:  however,  all  sensitivity  studies  will  be  carried  out  using  90  days  as  the  maximum 
number  of  days  the  plume  can  descend  while  retaining  its  integrity. 

B.  MODEL  PARAMETERS 
1.  Environmental  variables 

In  order  to  examine  the  environment's  influence  upon  possible  polyna-  related 
production  of  halocline  and  deep  waters,  certain  environmental  variables  were  individ¬ 
ually  varied  during  model  simulations.  The  sensitivity  of  the  model  to  a  variable's  al¬ 
teration  is  defined  by  the  dependence  of  the  plume's  depth  after  90  days  upon  the 
particular  variable,  relative  to  a  standard  case.  The  relationship  between  the  depth  of 
the  plume  and  a  single  environmental  variable  will  be  reported  in  Chapter  3.  The  vari¬ 
ables  to  be  investigated  are  listed  in  Table  I  on  page  14. 
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Figure  2.  Schematic  representation  of  the  coupled  model.  As  the  sheifwater  is  ex¬ 
posed  to  the  polyna,  a  high  density  plume  is  formed,  resulting  in  the 
plume's  descent  down  the  slope  until  it  is  neutrally  buoyant. 
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Table  1.  LIST  OF  ENVIRONMENTAL  INPUTS:  The  dependency  of  the  plume 
depth  after  90  days  on  each  of  the  variables  below  will  be  examined  in 
Section  3. 


Variable 

Definition. 

T.„ 

air  temperature 

r. 

wind  velocity 

s 

shelfwater  salinity 

r„ 

initial  shelfwater  velocity 

Sw 

shelf  width,  from  the  45  m  isobath 

initial  shelf  slope 

0, 

secondary  slope 

•V, 

Briint- Vaisala  frequency,  layer  1 

.Y, 

Briint-Vaisala  frequency,  layer  2 

AZ: 

thickness  of  layer  2 

2.  Environmental  conditions  for  the  standard  case 

To  achieve  consistency  with  pre-existing  works,  many  of  the  environmetal  con¬ 
ditions  for  the  standard  case  were  borrowed  from  Melling  and  Lewis  (1982),  Killworth 
and  Smith  (19S4).  and  Pease  (19S7).  The  standard  case  initial  values  are  listed  in 
Table  2  on  page  15. 

The  most  significant  departure  from  an  initial  value  given  by  previous  work,  is 
the  H,  value.  In  this  paper,  the  frazil  ice  collection  thickness  is  taken  to  be  .2  meters. 
Pease  (19S7)  shows  a  higher  collection  thickness  results  in  longer  polyna  opening  time 
and  larger  polyna  width  for  the  same  wind  speed  and  air  temperature.  Weeks  and 
Ackley  (CRREL  Monograph  S2-1.  1 9 S 2 )  give  higher  values  for  frazil  ice  collection 
thicknesses  in  a  wind  driven  scenario  vice  .01  to  .1  meters  in  quiescent  conditions.  In 
addition  to  the  frazil  ice  collection  thickness,  the  salinity  of  the  ice  is  set  to  7  ppt  (Cox 
and  Weeks,  1974)  instead  of  5  ppt  used  in  Killworth  and  Smith  (1984). 
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Table  2.  INITIAL  VALUES  USED  FOR  THE  STANDARD  CASE. 


Variable 

Value 

Definition 

T,„  "C 

-15. 

air  temperature 

T,.  "C 

-1.8 

water  temperature 

I sec1 

15. 

wind  velocity 

f  *.  m  sec-1 

3  °  o  r,„ 

ice  floe  velocity 

II .  m 

>  1 

frazil  collection  thickness 

o. 

5.67a*  10- 

*  Stephan- Boltzmann  constant 

e. 

0.95 

emissivity  of  the  air 

pckgnr2 

1.30 

air  density 

p?.  kgnr2 

1.026-vlO3 

seawater  density 

p.kgnr 3 

0.95xl03 

ice  density 

c 

2.0-vlO-3 

sensible  heat  coefficient 

C,  Jdeg-'kg-' 

1004 

specific  heat  of  air 

<X.  iv»r2 

301 

longwave  radiation  upward 

L.  Jkg-' 

3.34x10s 

latent  heat  of  fusion 

b, .  m 

10. 

flow  thickness 

u/,;  .  m 

1000. 

flow  width 

1  m  sec-1 

.04 

flow  speed 

/?  .  degrees 

29 

flow  direction 

/.  sec-1 

1.3Sxl0-* 

Coriolis  parameter 

.V.  sec-1 

0.0316 

Brunt- Vaisaia  frequency,  layer  1 

N 

0.0077 

Brunt- Vaisiila  frequency,  layer  2 

N 

0.001 

Brunt-Vaisala  frequency,  layer  3 

AZ,  .  meters 

195 

Layer  2  thickness 

0. 

0.5.rl  O'3 

initial  bottom  slope 

e3 

5xl0-3 

secondary  bottom  slope 

K 

0.01 

drag  coefficient 

E 

0.072  sin  6  sin  p 

entrainment  coefficient 

S  .  ppt 

32.5 

initial  shelfwater  salinity 

s,  ■  ppt 

7.0 

Salinity  of  frazil  ice 

15 


III.  RESULTS 


A.  STANDARD  CASE 

In  order  to  examine  the  environmental  influence  on  the  possible  polyna  production 
of  haloeline  and  deep  waters,  a  standard  case  was  defined  with  values  which  are  typical 
of  Arctic  Ocean  observations.  The  offshore  wind  velocity  requires  a  considerable 
meteorological  event  providing  winds  higher  than  the  monthly  mean  winds.  However, 
the  value  of  15  m  s  for  offshore  wind  velocity  is  within  observed  wind  speeds  for  Arctic 
storms  reported  in  previous  papers  (Schumacher,  et  al.  (19S3)  and  Aagaard.  et  al. 
(1985)).  The  standard  case  plume  obtains  a  salinity  of  34.98  ppt  after  transversing  the 
polyna.  which  corresponds  to  a  density  of  1032.93  gmjm 3  and  a  negative  buoyancy  of 
.0643  m/s 2  .  It  is  therefore  not  suprising  when  the  plume  penetrates  the  same  Arctic 
pycnocline  used  in  Melling  and  Lewis  (1982)  to  a  depth  of  436  meters.  Figure  3  on  page 
17  shows  the  velocity  of  the  plume  as  it  descends  the  slope.  Due  to  the  plume's  high 
negative  buoyancy,  the  plume  accelerates  slightly  on  the  initial  slope  then  rapidly  in¬ 
creases  in  speed  after  the  shelf  break  at  55  meters.  The  maximum  velocity  is  0.30  m  s 
reached  at  a  depth  of  60  meters.  The  mass  flux,  AV,  increased  by  over  1600%  empha¬ 
sizing  the  dependence  of  the  entrainment  rate  on  the  velocity  of  the  plume.  At  the  end 
of  90  days,  the  standard  case  plume  has  a  velocity  of  0.017  m  s  and  a  negative  buoyancy 
of  0.000482  ml  sec2.  The  plume  has  travelled  253  kilometers  along  the  shoreline  and  96 
kilometers  down  slope.  The  flow  path  of  the  plume  can  be  seen  in  Figure  4  on  page 
18.  Overall,  the  plume  in  the  standard  case  behaves  similarly,  although  with  different 
magnitudes,  as  the  plumes  in  Melling  and  Lewis  (1982).  With  this  in  mind,  sensitivity 
studies  of  the  impact  of  individual  environmental  parameters  on  the  90-day  plume  depth 
can  be  examined.  Sensitivity  experimental  results  typically  include  10-20  data  points  with 
increased  resolution  in  areas  of  abrubt  change. 

B.  SENSITIVITY  EXPERIMENTS 

The  polyna's  maximum  width  is  dependent  on  the  growth  rate  of  ice,  the  advection 
rate  of  ice,  and  the  collection  thickness  of  ice.  Air  temperature  and  the  wind  speed  af¬ 
fects  the  growth  rate  of  ice,  while  the  advection  rate  is  proportional  to  the  wind  speed. 
Lowering  the  air  temperature  results  in  a  monotonic  increase  in  the  plume's  depth  at  90 
days  as  shown  in  Figure  5  on  page  20.  This  is  expected  due  to  the  increased  production 
of  frazil  ice  resulting  in  higher  salinization  of  the  water  column.  Although  Pease  (19S7) 
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Figure  3.  The  standard  case  plume's  velocity  is  shown  with  depth:  Note  the  rapid 
increase  in  velocity  as  the  plume  crosses  the  shelibreak  at  the  55  meter 
isobath  to  a  steeper  slope. 

has  shown  the  polyna  size  decreases  with  lowering  air  temperature  at  a  given  wind  speed, 
the  plume's  salinity  (hence,  density)  increases  most  rapidly  in  the  first  10  kilometers  of 
downslope  travel.  Here  the  plume  will  be  subjected  to  significantly  higher  salinization 
due  to  higher  ice  growth  rate  and  lower  areal  flux  rate.  Although  the  polyna  size  is 
smaller  with  lowering  air  temperatures,  the  majority  of  the  salinization  has  already  oc¬ 
curred  in  the  first  10  kilometers.  Indeed,  the  polyna  shows  a  decreasing  salt  contribution 
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Figure  4.  The  standard  case  plume  non  path  in  the  x-y  coordinates:  The  effects 
of  Coriolis  can  be  seen  as  the  plume's  angle  P  decreases  during  the 
plume's  descent. 


IS 


after  the  first  10  kilometers  with  decreasing  air  temperature.  This  can  be  seen  in 
Table  3  on  page  19.  This  decreasing  salt  contribution  to  the  water  column  is  due  to  the 
polyna's  smaller  width. 


Table  3.  SALINITY  CONTRIBUTION  OF  POLYNA  WITH  DECREASING  AIR 
TEMPERATURE:  The  table  shows  the  initial  salinity  of  the  plume.  S,.„, 
after  the  first  10  kilometers  of  downslope  travel,  the  final  salinity  of  the 
plume.  S„  after  travelling  the  entire  width  of  the  plume.  A'.,  and  the  de¬ 
creasing  salt  contribution  of  the  polyna  during  the  distance  Xp  —  \0km  with 
decreasing  air  temperature.  T..r. - 


s„„ .  ppt 

5,  .  ppt 

AS  ,  ppt 

T  °  C 

1  cir* 

33.36 

34.36 

1.0 

-5.0 

33.75 

34.67 

0.92 

-10. 

34.14 

34.9S 

0.S4 

-15. 

34.54 

35.30 

0.76 

-20 

34.93 

35.61 

0.6S 

-25 

35.32 

35.92 

0.60 

-30 

35.72 

36.23 

0.51 

-35 

On  the  other  hand,  if  air  temperature  is  held  constant,  and  offshore  wind  speed  is 
increased,  the  polyna  grows  in  size  due  to  increased  ice  advection  rate.  The  initial  con¬ 
tribution  of  salt  during  the  first  10  kilometers  of  the  polyna  does  not  vary  however,  the 
increased  polyna  size  exposes  the  water  column  to  further  salinization.  This  results  in 
a  denser,  deeper  penetrating  plume  with  increasing  wind  speed.  (See  Figure  6  on  page 
21).  The  dependence  of  plume  depth  upon  the  wind  velocity  and  air  temperature  illus¬ 
trate  the  possible  atmospheric  influences  upon  the  Arctic  Ocean's  salinity  characteristics. 
The  frequency  and  severity  of  Arctic  storms  which  would  be  conducive  to  high  density 
plume  formation  is  a  consideration  in  the  calculation  of  production  rates,  but  cannot 
be  discussed  quatitatively  using  this  simplified  model. 

The  effect  of  varying  the  initial  salinity  of  the  shelfwater  on  plume  depth  is  shown 
in  Figure  7  on  page  22.  The  results  of  this  experiment  show  deeper  penetration  with 
higher  initial  salinity.  This  is  not  suprising  due  to  the  significant  role  of  salinity  upon 
water  density  near  the  freezing  point.  The  plume  depth  after  90  days  is  greater  than  the 
depths  reported  by  Melling  and  Lewis  (1982).  even  with  relatively  low  salinity  values. 

Initial  shelfwater  velocity  defines  two  parameters  within  the  model,  the  areal  flux 
used  in  computing  the  new  salinity  of  the  water  column  and  the  initial  velocity  of  the 
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Figure  5.  Effect  of  decreasing  air  temperature  on  the  plume's  90-day  depth:  With 
decreasing  air  temperature,  the  increased  ice  growth  rate  results  in  higher 
initial  salinization  of  the  water  column. 

plume.  The  higher  the  initial  velocity,  the  less  saline  the  plume.  However,  when  using 
the  plume  depth  at  90  days  as  a  standard  criteria  to  examine  the  results,  the  higher  ve¬ 
locity  plumes  reach  deeper  depths  at  90  days  titan  the  denser,  but  lower  velocity  plumes. 
If  a  speed  criteria  were  used,  then  there  would  be  a  negative  correlation  between  the 
plume  depth  and  initial  shelfwater  velocity.  Instead,  Figure  S  on  page  23  shows  a  pos¬ 
itive  correlation  between  plume  depth  and  initial  shelfwater  velocity. 
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polyna's  increased  width  with  increasing  wind  speed  results  in  a  higher 


density  plume  due  to  longer  exposure  of  the  water  column  to  the 
polyna's  brine  rejection.  This  results  in  a  deeper  penetrating  plume. 

Similarly,  the  parameters  of  shelf  slope  and  secondary  slope  appear  to  have  contra¬ 
dictors-  results.  Figure  9  on  page  24  shows  as  the  initial  shelf  slope  is  lessened,  the 
plume  penetrates  to  deeper  depths  after  90  days  of  travel.  In  this  experiment,  the  shelf 
width  was  fixed  at  20  kilometers  after  the  45  meter  isobath.  This  deeper  penetration 
appears  to  be  due  to  a  lower  entrainment  parameter,  which  is  a  function  of  the  shelf 
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Figure  7.  The  effect  of  initial  shelfwater  salinity  on  the  90-day  plume  depth:  With 
increasing  shelfwater  salinity  the  plume  attains  a  higher  density  due  to 
increased  salt  rejected  by  ice  growth  added  to  the  shelfwater. 

slope.  In  contrast.  Figure  10  on  page  25  shows  deeper  penetration  with  an  increased 
secondary  slope.  Although  the  entrainment  parameter  still  increases  with  slope,  the 
90-day  plume  depth  is  limited  by  the  slope  itself  Since  90  days  was  considered  a  rea¬ 
sonable  time  scale  for  the  plume  to  maintain  its  integrity,  deeper  plume  penetrations 
would  most  likely  occur  with  steeper  slopes  after  the  shelf  break,  i.e.  in  submarine 
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Figure  8.  Effect  of  varying  initial  shelfuater  velocity  on  90-day  plume  depth:  An 
increased  initial  velocity  results  in  a  deeper  penetrating  plume  at  90  days. 

canyons.  Varying  the  shelf  width  after  the  45  meter  isobath  has  little  effect  on  the 
90-day  plume  depth.  This  is  shown  in  Figure  1 1  on  page  26. 

Examining  the  stratification  of  the  Arctic  Ocean  is  accomplished  by  varying  the 
Brunt- Yiiisala  frequency  of  Layers  1  and  2.  The  deep  Layer  3  remains  constant  at  a 
value  used  in  Killworth  (1977)  for  deep  waters.  Layer  1  is  the  initial  pycnoc.line  and  re¬ 
presents  the  transition  from  the  fresh  surface  layer  to  the  saltier  halocline.  Figure  12 
on  page  27  shows  the  90  day  plume  depth  is  little  changed  even  with  a  substantially 
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Figure  9.  This  figure  shows  the  effect  of  the  initial  shelf  slope  on  the  90-day  plume 
depth:  Due  to  the  entrainment  parameter  being  a  function  of  slope,  a 
steeper  slope  results  in  higher  entrainment  of  ambient  water,  thus  more 
rapid  loss  in  negative  buoyancy. 

weakened  initial  stratification  with  A"2  =  2.0.tl0~4  sec2.  The  short  length  of  time  and 
space  which  the  plume  is  exposed  to  this  layer  is  probably  the  cause  of  its  diminished 
impact  on  the  90-day  plume  depth. 
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Figure  10.  The  effect  of  varying  the  secondary  slope  on  the  90-day  plume 
depth:  Even  with  a  higher  entrainment  rate,  the  plume  reaches  a 
deeper  depth  with  a  steeper  slope,  due  to  the  physical  boundary  of  the 
ocean  bottom. 

On  the  contrary.  Figure  13  on  page  28  shows  the  dramatic  difference  in  the  90-day 
plume  depth  with  small  changes  in  the  Brunt- Vaisala  frequency  in  this  2-15  meter  thick 
layer.  Additionally,  if  the  second  layer's  thickness  is  decreased,  the  plume  will  penetrate 
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Figure  11.  The  effect  of  varying  the  shelf  width  (after  the  45  m  isobath)  on  the 
90-day  plume  depth:  Varying  the  shelf  width  has  little  effect  on  the 
90-day  plume  depth. 

to  deeper  depths.  This  is  illustrated  in  Figure  14  on  page  29.  Note  the  discontinuity  in 
the  slope  of  the  curve  is  caused  when  the  plume  no  longer  pentrates  the  second  layer. 
Next,  the  standard  case  is  compared  in  detail  with  a  "submarine  canyon"  simulation. 
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Figure  12.  The  effect  of  varying  the  initial  pycnociine  strength  on  the  90-day  plume 
depth:  The  initial  pycnociine.  between  45  and  55  meter  depths,  has  a 
small  elTect  on  the  90-day  plume  depth  due  to  the  small  thickness  of 
this  layer. 

C.  THE  "SUBMARINE  CANYON"  CASE 

The  "submarine  canyon"  case  is  a  modification  of  the  standard  case  using  a  steeper 
secondary  slope,  lower  Layer  2  Briint-Vaisalii  frequency  value,  and  a  decreased  Layer  2 
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Figure  13.  The  effect  of  varying  Layer  2  Brunt- Yiiisala  frequency  on  the  90-day 
plume  depth:  Due  to  the  plume's  longer  exposure  time  in  this  layer, 
changes  in  Briint-Vaisala  frequency  tend  to  have  a  more  dramatic  effect 
on  the  vertical  motion  of  the  plume. 

thickness.  In  addition,  the  angle  /?.  is  forced  to  the  initial  value  of  29  degrees  to  simulate 
the  steering  of  the  plume  by  topography.  The  initial  values  used  in  the  "submarine 
canyon"  case  are  shown  in  Table  5  on  page  31.  Forcing  the  angle  /?  towards  its  initial 
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thickness,  the  90-day  plume  depth  becomes  constant. 


value  results  in  the  plume  travelling  down  the  V  axis  further  and  less  down  the  X  axis. 
Of  course,  the  value  chosen  for  forcing  /?  will  alter  the  X-Y  flow  path  of  the  plume.  The 
How  path  of  the  "submarine  canyon"  case  is  seen  in  Figure  15  on  page  32.  The  plume, 
descends  to  a  depth  of  over  1300  meters  in  90  days,  has  moved  207  kilometers  along¬ 
shore  and  160  kilometers  off-shore.  In  Figure  16  on  page  33.  the  "submarine  canyon" 
plume  is  seen  to  reach  a  maximum  velocity  of  0.36  m  s  whereas  the  standard  case  plume 
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has  a  maximum  velocity  of  0.30  m  s  (Figure  3  on  page  17).  The  effect  of  increased 
slope,  the  angle  /?,  and  velocity  on  entrainment  is  seen  when  comparing  the  volume  flux 
of  the  "submarine  canyon"  plume  (Figure  17  on  page  34)  and  the  standard  case  plume 
(Figure  18  on  page  35). 

D.  RESULTS  USING  OBSERVED  TEMPERATURE  AND  SALINITY  FIELDS 

The  values  for  this  experiment  are  the  same  as  the  standard  case  except  the  Briint- 
Vaisalii  frequency  is  computed  using  observed  values  of  temperature  and  salinity.  The 
temperature  and  salinity  records  are  from  Ostlund  (19S7)  and  Aagaard  et  al.  (1981) 
which  were  taken  in  the  Canadian  Basin.  The  results  of  the  experiment  are  shown  in 
Table  4.  The  90-day  plume  depth  and  maximum  velocities  are  very  similar  to  the  results 
from  the  Layer  2  Briint-Vaisala  frequency  sensitivity  experiment. 

In  an  attempt  to  model  a  current  as  observed  in  Aagaard  et  al.  (1985).  a  plume 
width  of  20  kilometers  was  used  together  with  a  steep  secondary5  slope  and  the 
temperature  salinity  fields  reported  in  Aagaard  et  al.  (1981).  The  plume's  velocity  with 
increasing  depth  is  shown  in  Figure  19  on  page  36.  Aagaard  et  al.  (1985)  reported  a  25 
kilometer  wide.  15  meter  deep  plume  with  a  mean  velocity  of  0.45  m  s.  The  simulation 
yields  a  velocity  of  0.30  m  s  before  the  shelf  break  and  0.92  m,  s  while  descending  the 
steep  secondary  slope. 


Table  4.  THIS  TABLE  COMPARES  SOME  RESULTS  FROM  USING  ACTUAL 
T-S  RECORDS:  The  AIWEX-2  and  AIWEX-8  stations  are  from 
Ostlund(  19S7)  and  the  third  record  is  taken  from  Aagaard  et  al.  (19S1). 


Variable 

AIWEX-2 

AIWEX-8 

Aagaard 

^  max  m  S 

.397 

.3S9 

.395 

90-day  depth,  m 

217 

213 

305 

X  distance,  km 

311 

312 

363 

'll'  distance,  km 

380 

375 

478 

30 


Table  5.  INITIAL  VALUES  USED  FOR  THE  "SUBMARINE  CANYON" 
CASE:  The  highlighted  values  are  those  parameters  which  differ  from  the 
standard  case. 


Variable 

Value 

Definition 

T, ...  C 

-15. 

air  temperature 

T„  .  "C 

-l.S 

water  temperature 

tn  sec*1 

15. 

wind  velocity 

V.  m  sec*1 

3  %  i;„ 

,  ice  tloe  velocity 

a .  m 

*> 

frazil  collection  thickness 

a,  lVnr:deg~4 

5.67x10-* 

Stephan- Boltzmann  constant 

0.95 

emissivity  of  the  air 

p7kgm~2 

1.30 

air  density 

kg  nr* 

1.026.vl03 

seawater  density 

p.kgnr* 

0.95.V  103 

ice  density 

C, 

2.0-rlO-3 

sensible  heat  coefficient 

C.  Jdcg-'kg-1 

1004 

specific  heat  of  air 

Q:u.  IVnr2 

301 

longwave  radiation  upward 

L.  Jkg 

3.34-vlO5 

latent  heat  of  fusion 

/?.,  .  m 

10. 

flow  thickness 

;r,  ,  m 

1000. 

How  width 

( m  sec*1 

.04 

How  speed 

ft  .  degrees 

29 

flow  direction 

/.  sec-1 

1.33.rlO-J 

Coriolis  parameter 

.V.  sec-1 

0.0316 

Brunt-Vaisala  frequency,  layer  1 

N 

0.006236 

Brunt-Vaisala  frequency,  layer  2 

A Z,  .  meters 

145 

Layer  2  thickness 

X 

0.001 

Brunt-Vaisala  frequency,  layer  3 

0, 

0.5.Y10-3 

initial  bottom  slope 

0, 

9y10-3 

secondary  bottom  slope 

K 

0.01 

drag  coefficient 

E 

0.072  sin  6  sin  (I 

entrainment  coefficient 

5  .  ppt 

32.5 

initial  shelfwater  salinity 

5,  .  ppt 

7.0 

Salinity  of  frazil  ice 

31 


forced  to  the  initial  />  value  to  simulate  the  topographic  steering  within 


a  submarine  canyon. 
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Figure  16.  The  velocity  of  a  "submarine  canyon"  plume  with  increasing  depth:  The 
plume  has  the  characteristic  "jump"  in  velocity  after  the  shelf  break 
where  it  reaches  a  maximum  velocity  of  0.36  m.  s. 
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water  reaching  the  deeper  depths  is  shelf  water. 
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PLUME  DEPTH  ( M ) 


Figure  19.  Velocity  shown  with  depth  of  a  20  kilometer  wide  plume:  This  figure 
shows  a  wider  plume  can  penetrate  deeper  than  the  narrower  1  kilo¬ 
meter  plumes  used  in  the  previous  experiments. 
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IV.  DISCUSSION 


This  paper  has  attempted  to  provide  a  mechanism  in  which  high  salinity,  low  tem¬ 
perature  water  is  produced  and  transported  otT  the  continental  shelves  into  the  Arctic 
halocline  or  deeper  depths.  The  environmental  parameters  which  affect  the  production 
and  transport  of  this  cold,  salty  water  were  examined,  resulting  in  the  conclusion  that  a 
polyna  can  produce  water  with  a  salinity  greater  than  34.5  ppt  near  the  freezing  point 
which  can  be  transported  via  gravity  flow  plumes.  The  results  show  deeper  penetrations 
to  depths  greater  than  700  meters  is  possible  depending  on  the  set  of  environmental  pa¬ 
rameters.  This  indicates  polyna-produced  high  density  plumes  may  also  be  responsible 
for  the  ventilation  of  the  deep  Arctic  Basins,  specifically  the  Canadian  Basin. 

Ostlund  et  al.  { 1 9S7)  placed  the  residence  time  of  Canadian  Basin  Deep  Water  be¬ 
tween  500  and  800  years.  Using  the  mixing  ratio  of  2:1  (shelfwater  to  intermediate  wa¬ 
ter)  from  Aagaard  et  al.  (19S5),  shelfwater  with  a  salinity  of  35.1  ppt  requires  an 
off-shelf  flux  of  approximately  .0063  Sv  to  ventilate  the  Canadian  Basin  in  800  years. 
The  results  from  the  20  kilometer  wide  plume  experiment  give  an  annualized  flux  of 
0.0035  Sv.  In  other  words,  an  average  of  about  two  such  events  per  year  where 
polyna-produced  plumes  penetrate  below  the  temperature  maximum  over  the  past  800 
years  could  account  for  the  ventilation  of  the  Canadian  Basin.  Furthermore,  Ostlund 
et  al.  (1987)  states  that  the  deep  water  probably  consists  of  much  less  than  10% 
shelfwater  and  would  have  to  be  a  salinity  of  36.15  ppt.  The  "submarine  canyon"  case 
shows  the  plume  reaching  a  depth  of  over  1300  meters  consisting  of  only  1.3% 
shelfwater  at  this  depth.  Table  3  on  page  19  shows  under  extremely  cold  conditions,  a 
salinity  value  of  36.23  ppt  can  be  achieved  with  this  model.  Therefore,  it  is  reasonable 
to  expect  that  polyna-produced  saline  water  may  significantly  contribute  to  the  venti¬ 
lation  of  the  Canadian  Basin. 

Climatological  variations  which  would  affect  the  number  of  deep  penetrating  plumes 
are  worthy  of  discussion.  Besides  the  major  periods  of  global  ice  ages  and  their 
counterpart  warm  eras,  there  have  been  recorded  several  shorter  cold  and  warm  epochs. 
These  shorter  periods  coincide  with  weak  and  strong  wind  circulation,  respectively.  In 
the  20th  century,  both  a  warm  period  lasting  about  20  years  and  a  colder  period  begin¬ 
ning  around  1950  have  been  recorded.  In  each  of  these  two  eras,  there  has  been  a  sud¬ 
den  outflow  of  surface  layers  from  the  Arctic  Ocean.  Increased  ice  flow  was  recorded 
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in  1938  around  Iceland  during  the  warm  period,  while  an  increased  flow  of  cold,  low 
salinity  water  reached  Iceland  from  the  Arctic  in  the  1%0's.  Removal  of  the  colder, 
fresher  surface  layer  from  the  Arctic  Basin  would  reduce  the  vertical  temperature  and 
salinity  gradients  resulting  in  a  less  stable  vertical  structure  for  the  ocean.  In  effect,  an 
environment  similar  to  reducing  the  model's  Layer  2  Brunt-Viiisala  frequency  and  re¬ 
ducing  the  Layer  2  thickness  may  occur,  enhancing  deep  plume  penetrations.  In  addi¬ 
tion  to  possible  weakening  of  the  vertical  structure,  removal  of  the  ice  cover  leaves  large 
areas  of  open  water  where  heat  is  more  readily  removed  from  the  ocean.  In  these  large 
areas  of  open  water,  frazil  ice  grows  rapidly  and  abundantly  to  increase  the  brine  rejected 
into  the  water  column.  In  the  warm  period  of  the  early  20th  century,  the  extent  of  the 
Arctic's  semi-permanent  ice  pack  was  greatly  reduced  (up  to  10%)  between  1920  and 
1938.  The  reduction  in  the  semi-permanent  ice  pack  is  attributed  to  higher  ice  removal 
rates  (increased  advection)  due  to  the  strong  wind  circulation  of  the  period.  In  Ostlund 
et  al.  (1987).  it  is  noted  that  the  deepest  waters  in  the  Canadian  Basin  may  be  the  re¬ 
mains  of  the  warm  climate  which  occurred  around  1000  A.D.  vice  any  cold  epochs  such 
as  the  "Little  Ice  Age”  in  the  late  17th  and  early  18th  centuries. 

In  summary,  periods  when  Arctic  wind  circulation  is  strengthened  resulting  in  open 
water  due  to  rapid  ice  divergence  coupled  with  a  weaker  vertical  Arctic  Ocean  stability 
seems  to  be  the  most  conducive  environmental  scenario  for  deep  plume  penetrations  to 
provide  the  additional  salt  source  in  the  Canadian  Basin. 

In  Killworth  and  Smith  (19S4),  the  heat  diffusivity  constant,  k  ,  required  to  achieve 
steady  state  balances  requires  some  peculiar  changes  in  magnitude  with  increasing  depth. 
A  value  of  k^3  x  10-4m\s_1  is  established  to  achieve  a  steady  balance  between  the  non- 
turbulent  shelfwater  plume  used  in  Killworth  and  Smith  (1984)  and  the  region  where  the 
plume  interleaves  (note:  the  density  of  this  plume  is  based  on  the  brine  rejection  due  to 
a  single  year's  ice  growth  over  a  continental  shelf  of  50  meters  average  depth).  To 
properly  supply  the  East  Greenland  current  outflow,  k  was  required  to  be 
so  x  10_-‘m2  sec-1.  And  finally,  to  balance  the  inflow  of  warm  Atlantic  water  and  upward 
advection  of  heat,  k  at  this  level  was  required  to  be  ^1.2  x  10'4m2  sec*1  If  the  normalized 
mass  flux  of  the  turbulent  plume  model  was  to  be  used  as  a  gauge  to  understand  the 
turbulent  diffusion,  this  thesis  would  present  some  interesting  results.  In  Melling  and 
Lewis  (1982),  the  normalized  mass  flux  was  of  order  ^10°.  In  the  standard  case,  the 
normalized  mass  flux  was  of  order  ^10l.  The  "submarine  canyon"  case  had  a  mass  flux 
of  order  ^102.  Could  it  be  that  the  three  above  cases  provide  a  insight  into  the  missing 
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physics  so  stated  in  Kilhvorth  and  Smith  ( 1 9S4 >?  Some  other  missing  physics  may  be  the 
double  diffusion  process  as  outlined  in  Carmack  and  Aagaard  (1973)  and  the  breaking 
of  internal  waves  (Perkin  and  Lewis.  1978).  Aagaard  et  al.  (1985)  suggests  the  incorpo¬ 
ration  of  turbulent  plume  dynamics  into  the  "filling  box"  model  of  Kilhvorth  and  Smith 
(19S4).  A  mechanism  providing  a  specific  low  temperature,  high  salinity  water  source 
has  been  developed  in  this  paper.  Although  the  mechanism  is  episodic,  as  suggested 
Aagaard  et  al.  (19S5),  (Kilhvorth  and  Smith  (1984)  had  earlier  used  the  term 
"spasmodic"),  due  to  the  dependence  on  environmental  variables,  flux  rates  of  this  highly 
saline  water  can  be  estimated  from  renewal  times  of  the  Canadian  Basin.  The  problem 
of  incorporating  the  turbulent  plume  in  the  "filling  box"  model  may  be  tackled,  but  is 
beyond  the  scope  of  this  paper. 

The  maintenance  of  the  Arctic  halocline  by  the  mechanism  presented  in  this  study 
alone  would  require  many  times  the  number  of  active  polynas  producing  saline 
shelfwater.  Using  the  3:2  mixing  ratio  proposed  by  Aagaard  et  al.  (1981).  a  minimum 
of  2.5  Sv  production  rate  of  cold,  saline  shelfwater  was  determined  necessary  to  replenish 
the  halocline  over  a  period  of  10  years.  With  the  flux  rate  of  the  20  km  wide  plume  (.0035 
Sv),  the  individual  mechanism  modelled  in  this  paper  would  require  200-300  coastal 
polynas  producing  2-3  high  salinity,  low  temperature  plumes  per  year  to  maintain  the 
Arctic  halocline.  Although  200-300  such  polynas  may  be  thought  prohibitive,  areas  of 
open  water  which  exist  during  the  onset  of  freezing  temperatures  in  the  late  fall,  as  well 
as  leads  which  continue  through  winter  over  as  much  as  10%  of  the  ice  pack,  could 
provide  a  rapid  source  of  salt  to  help  overcome  the  summertime  salinity  gradient.  In 
times  of  strong  vertical  structure,  these  polyna-produced  plumes  can  transport  low- 
temperature.  high  salinity  water  into  the  Arctic  halocline.  Overall,  the  polyna  seems  to 
be  an  important  source  of  brine  which  may  help  maintain  the  Arctic  halocline  and  ven¬ 
tilate  the  deep  basins. 
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